On the theory of MHD waves in a shear flow of a magnetized turbulent plasma 
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The set of equations for magnetohydrodynamic (MHD) waves in a shear flow is consecutively 
derived. The proposed scenario involves the presence of a self-sustained turbulence and magnetic 
field. In the framework of Langevin-Burgers approach the influence of the turbulence is described by 
an additional external random force in the MHD system. Kinetic equation for the spectral density 
of the slow magnetosonic (Alfvenic) mode is derived in the short wavelength (WKB) approximation. 
The results show a pressing need for conduction of numerical Monte Carlo (MC) simulations with a 
random driver to take into account the influence of the long wavelength modes and to give a more 
precise analytical assessment of the short ones. Realistic MC calculations for the heating rate and 
shear stress tensor should give an answer to the perplexing problem for the missing viscosity in 
accretion disks and reveal why the quasars are the most powerful sources of light in the universe. 
The planned MC calculations can be incorporated in global models for accretion disks and also in 
all other physical conditions where there is a shear flow in a magnetized turbulent plasma. It is 
supposed that the heating mechanism by Alfven waves absorption is common for many kinds of 
space plasmas from solar corona to active galactic nuclei and the solution of these longstanding 
puzzles deserves active interdisciplinary research. The work is illustrated by numerical calculations 
and by exact solutions for the time dependence of the magnetic field given by the Heun function. 

PACS numbers: 52.35.Bj, 95.30.Qd, 98.62. Mw, 97.10.Gz, optional 47.20.Ft 



I. INTRODUCTION 

The purpose of the present work is to give a de- 
tailed derivation of the stochastic magnetohydrodynamic 
(MHD) set of equations for a shear flow in a magnetized 
turbulent plasma. The meticulously performed analysis 
of the MHD system sets the basis for further Monte Carlo 
(MC) calculations devoted to reveal the basic physical 
phenomena in accretion disk plasmas: their heating and 
the origin of a large effective viscosity significantly ex- 
ceeding the bare plasma one. The paper is organized as 
follows: short historical remarks introducing to the read- 
ers the motivation for the current work are presented in 
Sec. mi Basic wave kinematics in a shear flow with a tran- 
sition from Eulerian to Lagrangian coordinates is given in 
Sec. IIIII Derivation and linearization of the full stochas- 
tic MHD set of equations with external random forcing 
by turbulence is described in Sec. IIVI Plane waves an- 
zatz in a shear flow is included for separation of variables 
and subsequent reduction from partial to ordinary differ- 
ential equations is performed in Sec. [V] The illustrative 
2D case is analyzed in Sec. [VI] Secular equation for the 
Alfven waves amplitudes is solved and the correspond- 
ing damping rate is obtained. Auxiliary problem for the 
period averaged energy of an effective oscillator under a 
white noise is considered. Short wavelength (WKB) ap- 
proximation is applied to the Alfven waves amplitude in 
Sec. IVIII and is it shown that the Alfven spectral den- 
sity obeys an effective Boltzmann equation. Sec. IVIIII 
treats the Langevin-Burgers MHD, modeling the influ- 
ence of the turbulence as a random external force in the 
momentum equation. Further on, speculations on the 
origin of a large effective viscosity are proposed. Conclu- 



sive remarks and future perspectives are briefly discussed 
in Sec. IIXI It is debated on the kind of numerical anal- 
ysis which has to be done in order to reveal the origin 
of the huge effective viscosity, observed in the accreting 
magnetized turbulent plasma. 

All the analytical calculations utilized in the current 
work, as well as some general concepts for heating, re- 
lated to the stochasticity of the investigated turbulent 
system, are laid out in five separate appendices. Matrix 
presentation for Lagrange-Euler transformations in the 
presence of a shear flow is reproduced in appendix [X] 
Linearization of the dynamic equations is performed in 
Appendix [B] Detailed derivation of the complete MHD 
set of equations for a magnetized plasma in a shear flow 
under the influence of a random noise is provided in ap- 
pendix [UJ Test examples of MHD waves with a restricted 
wave-vector orientation K y — 0, short wavelengths and 
small attenuation are presented. The shear rate depen- 
dent damping of Alfven waves propagating along the 
magnetic field lines is calculated in appendix [D] Several 
illustrative examples for the power rate due to stochastic 
heating of a Brownian particle, oscillator under a white 
noise and a free particle are applied in[E]to support the 
general consideration for a white-noise driven heating, 
introduced in sec. IVIIH 



II. MOTIVATION 

Matter accreting onto a compact object redistributes 
radially its angular momentum, dissipates gravitational 
energy and forms a powerful source of radiation through 
physical processes that still cannot be regarded as un- 
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derstood. Accreting plasma's temperature and pressure 
predict a molecular viscosity that is orders of magnitude 
too small to account for the observed radiation inten- 
sity. Thus the origin of a large effective viscosity in the 
magnetized turbulent plasma of the accretion disks is a 
significant and yet unsolved problem. We know how a 
laser emits light, how the luminescent lamp works, we 
know how flash a fire-fly in the summer nights, but we 
do not know yet what is the mechanism of glowing of the 
most luminous sources of light in the universe - we plan 
to reveal this long-standing problem. 

Frictional release of angular momentum is believed to 
be an important element in stars formation as well. In 
our solar system, for instance, the large mass concen- 
trated into the sun carries only two percents of the an- 
gular momentum, 1 while all the rest is associated with 
the planetary motion. Next to the Galileo assertion "Ep- 
pur si muove " (and yet it does move) , for half a century 
we still face the question "why does the Sun not move 
(rotate) faster". Presumably, at an early stage of the 
system's formation the proto-planetary disk has acted 
as a brake which slowed down the rotation of most of 
the disk's mass, so it could cluster and "ignite" as a star. 
Broadly, such a mechanism of angular momentum release 
appears responsible for the formation of compact astro- 
physical objects and possibly determines the universe in 
the way that we observe it today. In this context, there 
still exists a central open question of what physical pro- 
cess produce the friction forces that facilitate the stars 
formation. 

The purpose of the present work is to demonstrate a 
strong energy dissipation in a shear flow of a magnetized 
plasma based on a simple model. Despite the process's 
interdisciplinary nature, we will pursue a description that 
is built on first-principles physics. Our leading hypothe- 
sis is that in a shear flow of the almost inviscid plasma 
the Alfven waves get amplified^ and later have their en- 
ergy thermalized, i.e. "lasing" of Alfven waves is the 
basis of dissipation in accretion disks. We select an ap- 
proach which is statistical rather than fluid-mechanical, 
by employing a kinetic equation for the spectral density 
(proportional to the square of their amplitude) of the 
Alfven waves. In the spirit of the quantum mechanics 
the square of the amplitude is proportional to the num- 
ber of particles, which we will call "alfvons". Here the 
kinetic equation describes the dynamics of the alfvons 
population. The lasing of the media is analogous to the 
dynamics of ecological systems where a fast population 
growth is followed by a resettlement and high death-rate. 
Similarly here, the energy transferred from the shear flow 
first increases the alfvons density and is later dissipated, 
thus effectively raising the plasma's viscosity and resis- 
tivity. Our goal is to deduce a kinetic equation for the 
spectral density of the Alfven waves, as well as its solu- 
tion, and analyze the physical factors driving the process. 
Standard methods of quantum mechanics (quasi-classical 
approximation for short wave-lengths and perturbation 
theory to account for the small initial viscosity of the 



plasma) will be considered as adequate. As an initial 
step, we employ a simplified treatment of the turbulence 
in accretion disks - the Burgers approach and show that 
the turbulence's role in a model approximation may be 
reduced to the source term in the kinetic equation. We 
also assume that it is the turbulence that triggers the 
Alfven waves, which are later amplified by the shear flow. 
There are many examples where waves can be amplified 
but often conditionally one can mention lasers and in 
general case lasing processes in some unstable medium. 

The considered model problem involves convective in- 
stability and turbulence in the heated disk's plasma and 
describes a self-sustained scenario for heating in accre- 
tion disks. The velocity fluctuations serve as a random 
force that initiates the Alfven waves, which are then sub- 
ject to a large amplification. The energy carried by these 
alfvons is then transformed into heat through the molec- 
ular viscosity and Ohmic resistivity, thus also creating 
friction forces. As the heat generation is more intense in 
the middle of the disk's thickness, the consequent tem- 
perature difference drives the convective instability and 
the turbulent convection. Thus the process sustains itself 
and the gravitational energy transforms into heat. 

III. WAVE KINEMATICS IN SHEAR FLOWS 

To describe locally the motion of the accreting fluid, 
we choose the z axis along the velocity and the x axis in 
direction of the velocity's gradient. Thus for the back- 
ground shear velocity field we have 

= Ax, t/W = = 0. (1) 

The superscript (0) refers to an equilibrium laminar shear 
flow whose perturbations will be studied. 

A selected small element of the fluid is carried by the 
flow so that its z coordinate is a linear function of time 

z at (t) = Axt + z at (0). (2) 

Under such drift the coordinates x, y remain unchanged, 
i.e. x a t(t) = const, y a t(t) = const. For a hydrodynam- 
ical description we will use flow-following (Lagrangian) 
coordinates 

x = x, y = y, z = z — Axt. (3) 

These tilde variables f are the Cartesian coordinates 
of the "tagged" atoms at the initial moment z at (0) = 
z a t (t) — Axt. They express the initial position of the fluid 
element. The change of the variables ([3]), however, does 
not alter the projection in the xy plane. As the initial 
position of the "tagged" atoms is fixed we can consider 
the tilde coordinates as "frozen" in the fluid. 

Let us now consider a plane wave in the shear flow with 
an amplitude cx exp(ik • r). The requirement of a phase 
invariance 

k f = k r (4) 
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sets the transformation law 



(5) 



into 



which can be validated by a substitution of Eq. 
Eq. g). 

As the initial tilde coordinates are related to the frozen 
initial position of the atoms, the wave vector in tilde 
coordinates is also "frozen", k = const. This determines 
the evolution of the wave- vector in Cartesian coordinates 



k x (t) 



Atk 7 



h h 



(6) 



This time dependence of the wave vector k(f ) has purely 
kinematic origin and it is not related to the dispersion 
of the waves. Such a phenomenon is well-known in 
the acoustics of moving media, but it can affect even 
non-propagating spatial structures with "tagged" atoms. 
Even in this static case with exactly zero frequency the 
general formula for the wave-vector evolution Eq. ^ is 
applicable. 

In the next section we will apply these kinematic re- 
lations to the Alfven waves. We consider the Alfven ve- 
locity Va, defined by an external magnetic field Bo, the 
density of the fluid p and magnetic pressure p B 



Pe 



2po ~ 2 pVa > 



(7) 



and the shear parameter A with dimension of a frequency. 
In Gaussian system the magnetic permeability of the vac- 
uum is po — 47r; in SI for simplicity po = 4ttx 10~ 7 . For a 
kinematic description it is convenient to introduce dimen- 
sionless wave vectors K = (Va/A) k, K = (Va/A) k, as 
well as dimensionless time r = At; La = Va/A is the 
unit for length. The connection between the Eulerian 
components of the wave-vector and the components in 
the tilde Lagrangian system Eq. (jJJJ) reads as 



K x (t)=K x {0)-tK z , 
K v (t) = K y (0), K z (t) = K z (0). 



(8) 



With an appropriate choice of the initial time we can 
further set K x (0) = 0. Here we have used that the wave- 
vector in the Lagrangian (tilde in our notations) coordi- 
nate system is constant, K(r) = const. Lagrangian co- 
ordinates are also known as flow-following coordinates, 
which in this particular flow preserve the wave vector 
(similarly, a plane of "tagged" atoms is moved but not 
deformed). For matrix presentation of the considered re- 
lations see Appendix |A"1 



IV. LINEARIZED STOCHASTIC MHD 

Consider the laminar component of the velocity field 
V(t, r) of an incompressible flow of an accreting plasma 
with a constant density p = const. In the presence of an 
external magnetic field B(t, r) the velocity V(t, r) evolves 



according to the Navier-Stokes equation with a Lorentz 
force 



p(d t + V- V)V = -Vp + rjAV + j x B 



(9) 



where p is the pressure and rj is the viscosity^. The term 
P(t, r) phenomenologically describes the Reynolds and 
Maxwell stresses associated with the turbulence. We will 
return to this term later, when we analyze the evolution 
of perturbations in the shear flow. 
The current j is given by Ohm's law 



j=a(E + VxB), 



(10) 



where a = 1/g is the electrical conductivity and (E + 
Vx B) expresses the effective electric field acting on the 
fluid; the velocity V <C c is nonrelativistic. 

For low frequencies j ^> £o|9*E|, so it is possible to use 
the magnetostatic approximation 

rotB = M oj. (11) 

We substitute here the current from Eq. (fT0|) and obtain 

E = -V x B + v m rot B, (12) 

where v m = £oc 2 g is the magnetic viscosity. In Gaussian 
system e — If Am while in SI e = 1/(4tt x 10~ 7 c 2 ). 

Here we will remind some basic properties of the 
plasma. We suppose that the frequency of ion-ion colli- 
sions t>a is much bigger than the ion cyclotron frequency 



v-n > wb 4 , oJBi = eB/M, 



(13) 



where M is the mass of the ion, see Ref. [6|, sees.: 41- 
43, 58. For weak fields we can neglect the influence of 
the magnetic field over the viscosity and resistivity; see 
Ref. @|, eqs. (43.8-43.10) 



OAT, 



5/2 



V 



a 



0.6T, 



3/2 



e 4 7V p M 1 /2£ p ' e 2 mV2£ e 



0.6 * 4tt T t 



3/2 ' 



e = 9 e /47T£ , 



(14) 
(15) 



where m is the electron mass, q e is the electron charge, 
N is the number of electrons per unit volume, C is the 
Coulomb logarithm and ro is the Debye length 



C 



I'D 



T 



4ire 2 N 



For high enough temperatures T 3> Tk m , where 

T£ m = 4nmc 2 e 6 N£\/^-, 
V m 

the kinematic viscosity dominates v m . 



(16) 



(17) 
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The evolution of the magnetic field is governed by the 
other Maxwell equation rotE = — d t B, so that the in- 
duction equation reads 

d t B = rot (V x B) - u m rot rot B. (18) 

For an incompressible flow div V=0, taking into account 
that divB=0, the equation above takes the form 

(d t + V-V)B = (B-V)V + ^ m AB, (19) 

where we have used the general relation 

rot(V x B) = VdivB + (B • V)V - B divV - (V • V)B. 

We will use the dynamics equations Eqs. ([HI [T9|) to ana- 
lyze the propagation of MHD waves in a shear flow. Let 
the velocity V be a sum of equilibrium shear velocity 
Tj(°)(r) and a small perturbation u(t,r), for which we 
are going to derive linearized wave equations 



V = U (0) + u. 



(20) 



The same representation we suppose for the magnetic 
field and the pressure 



B = B +B, p = Pa +p . 



(21) 

As mentioned in section IIII| the z-axis of the coordi- 
nate system we choose along the shear velocity and x-axis 
along the velocity gradient 

U (0) = (0,0, Ax) = Axe z , e 2 = (0,0,1). (22) 

For this choice of the coordinates the vorticity 

rotU (0) = (0,-A0) = -Aety, e v = (0,1,0), (23) 

is along y-axis. 

The differential rotation in the accretion disks stretches 
the frozen-in magnetic field lines and eliminates the cross- 
flow magnetic field. In the present work we will consider 
the magnetic field in the plasma to be parallel to the 
shear flow in z-direction 

B = (0,0, S ) = B e z . (24) 

We will consider a general perturbation case 

u(t,r) = (u x ,u y ,u z ), B'(t,r) = {B' x ,B' y ,B' z ). (25) 

Neglecting the small quadratic terms (u • V)u and (u • 
V)B in the substantial derivatives in Eqs. ([9J [T9| the 
linearized evolution equations read (see Appendix |B|) 

Vp f 

(dt + Axd z )u = Au x e, -i 

P P 



B 



'd z B x - d x B z 



jh ^d z B' y -d y B' z I Au. (2(i) 

(d t + Axd z )B = B d z u 

+AB' x e z + v m AB' . (27) 



For the density of the random force that models the tur- 
bulence we assume a white noise correlator 

(F(t l5 n)F(< 2 , r 2 )) = f p 2 6(t 1 - t 2 )<y(n - r 2 )l, (28) 

parameterized by the Burgers parameter T. In such a 
way we derive in Eulerian coordinates r = (x, y, z) a sys- 
tem of partial differential equations. The transition to 
Lagrangian f = (x,y, z) variables Eq. ([3]), however, re- 
duces it to a system of ordinary differential equations. 
To Eq. ((3]) we add also t = t and perform the change of 
variables 

u(t,r) = u(f,f), B'(i,r) =B(t,r), F(t, r) = F(t, f), 

(29) 

supposing that (t, r) are 4 independent variables as (t, r) 
are. From Eq. ([3]) one can easily derive the rules for the 
change of variables in the derivatives 

d t = df-Axds, d x = d x ~Atd z , d y = d y , d z = d~ z . (30) 

The main advantage of this transition to Lagrangian vari- 
ables is that the substantial time derivative dt + Axd z = 
di in the MHD set of equations Eq. (|26j) no longer de- 
pends on the spatial coordinates. 



V. SEPARATION OF VARIABLES 

For a system with constant (space independent) coef- 
ficients the solutions of the equations are plane waves. 
That is why using the phase invariance Eq. Q one can 
seek a solution in the form 

u = u^(t) exp(ik • f) = Uj c (t) exp(ik(t) ■ r), (31) 
B' = B~ (?) exp(ik ■ f) = B'^(t) exp(ik(t) ■ r), (32) 

where Uy.(i) = u^(t) and B^(t) = B~(£) are the time 
dependent amplitudes of the plane waves. For now we 
suppose that the variables are complex- valued. 

For clarity we will perform the change of variables 
gradually and will analyze the result for each term in- 
dividually. First we will change the variables only in the 
substantial time derivatives so that the system Eq. ([26]) 
takes the form 

F 

dfU + Au x e z = 

P 

Vp Bn ( dzB T ~ d * B f\ 
-— + — d z B y - d y B z + ^Au, (33) 
P PoP\ Q I 

5 t -B' =B d z u + AB x e z + u m AB', (34) 
d x B x + d v B y + d z B z = 0, (35) 
d x u x + d y u y + d z u z = 0. (36) 



To emphasize the physical meaning we have used La- 
grangian variables, however one can start with the simple 
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relation which both automatizes and ensures the separa- 
tion of variables 

(d t + Axd z ) (u(t) exp |z (k x — Atk z )x + k y y + k z z 

= exp|i (k x — Atk z )x + k y y + k z z j d t u(<). (37) 

For the plane-wave in Eq. (f3"Tj) the substantial derivatives 
A = d t + V ■ V (38) 
after linearization take the form 

D t [u £ (t)exp(tk(t)T)] 
« exp(ik(i) ■ r) [d t -u fc (i) + ^ k (t)e z ] , (39) 
A [Bj(t) exp(zk(i) • r)] » exp(ik(t) • r)dfBj(t). (40) 

As mentioned earlier, in Eqs. (|391 HD|) we have neglected 
the quadratic terms, which describe a wave-wave inter- 
action. 

Now we can substitute here the right-hand side of 
the supposed plane form of the waves from Eq. (|31[). 
where they are presented in Eulerian coordinates. The 
latter are more convenient to calculate the gradients 
V = ik(t) and Laplacians A = —k 2 (t) but we have 
to remember that according to Eq. ([6]) the wave-vectors 
k(i) = (t4/Va)K(t — At) are time dependent 

k(i) = {k x -Atk z , k y , k g ), k 2 (t) = k 2 x + k 2 y + k 2 z . (41) 

Additionally, to obtain a system of ordinary differential 
equations for the amplitudes of the MHD waves, one can 
eliminate the pressure as it is done in Appendix [Cl We 
suppose plane waves for the velocity and the magnetic 
field 



tcm of equations: 



u(t, r) = u^(t) exp(ik(t) • r), 
B'(t,r)=Bj(t)exp(ik(t).r), 



(42) 
(43) 



where we have introduced several notations that will be 
used later 

u kW = -*Uk(*) = -iVAV^r), (44) 
B k (t) = S b £ (r). (45) 

For the sake of brevity the Lagrange wave- vector indices 
k will be further suppressed; for example, taking the real 
part from Eqs. (J42j |t3|) we have 



u(t,r)=V^t7(r)Bin(k(t).r), 
B'(i,r) = B b(r)cos(k(t) -r). 



(46) 
(47) 



For the real dimensionless amplitudes v and b after 
some algebra given in the Appendix [C] we derive the sys- 



d T v x 

d T v y 

d T b x 
d T b y 

v. 



2K K 

-j^v x -K z b x -v' k K 2 v x + gf x , (48) 
2K K 

— ^V X - K Z by ~ V ' k K 2 Vy + hfy, (49) 

K z v x -v' m K 2 {r)b x , (50) 

K z v y ~ V ' m K 2 (r)by, (51) 

(53) 



where 



9\J 



K v + K z JK 2 (t)+K 2 

, h(r) = V xK \ z , (54) 



K{t) ' v ' K{t) 
K x = —tK z , K y — const, K z = const, (55) 

^/k 2 + (1 + t 2 )K 2 , (56) 



K(r) = ^K 2 + K 2 + K 2 



A 



A 



vf k ' Vm ~ V 2 



A 



V 2 AV 



(57) 



v = v k + v m , v' = v' k + v' m , (58) 
(/x(n)/ x (r 2 )) = (f y (n)f v {n)) = r5(n ~ (59) 

and V is the volume of the system. 

In the simplest case of zero shear A = and dissipa- 
tion v k — v m — we have a system with constant coef- 
ficients which gives in dimensionless form the dispersion 
of Alfven waves [f|, sec. 69, problem 



V a i 
cu A (k) = -i|k.B |-- l /k 2 , V A 

£>0 A 



By 
Bo 



v gr (k) = ^i = y AS ig n (k.B 0/B 



(60) 



(61) 



The wave- wave interaction is small and negligible only if 
the dimensionless wave components of the velocity and 
magnetic fields are sufficiently small v 2 , d 2 < 1. 



VI. TWO DIMENSIONAL WAVES 

To investigate the influence of the shear flow on the 
evolution of the MHD waves, here we will concentrate on 
the two-dimensional case k y — which gives a complete 
separation of the variables. For this special case K y = 
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the dynamic equations take the form 



d T v x 


= a(r)v x - K z b x 






-v' k K 2 {T)v x +g( T )f x {T), 


(62) 


d T v y 


= -K z b y - V ' k K 2 (r)v y + f y (r), 


(63) 


A h 

& T O x 


— J^zV x ~ V m l\ (T)O x , 




d T b y 


= K Z Vy - v' m K 2 {T)by, 


(65) 


V z 


K x 

K z Vx ' 


(66) 


b z 


K x , 


(67) 



where 



„ KJt)K z 
a(r) ee 2- 



2r 



KI(t) + KI 



K 2 (r) = (l+r 2 )Kl g(r) = 



(68) 
(69) 



We will start our analysis with the case of short wave- 
lengths in a dissipationless and a fluctuation-free regime. 
For a highly conducting plasma v m w with a negligible 
kinematic viscosity v k ss (conditionally, we may say su- 
pcrfluid and superconducting plasma) the system above 
yields 



d T v x = a(r)v x - K z b x , 
d T b x = K z v x . 



(70) 
(71) 



We differentiate Eq. {7DJ) with respect to time, neglect 
(for \K Z \ 1) the d T a term, and substitute d T b x from 
Eq. (I7ip which implies 



v — cx(t)v + K 2 v fa 0, 



(72) 



where v = v x and the dot stands for the dimensionless 
time t derivative. In the original variables this equation 
reads 

d K,k + ls(t)dtu x ^ + u\u x £ « 0, (73) 
where the frequency of the Alfven waves 

w A (k) = V A \k z \ = V A k\ cos^|, (74) 
depends on the angle 

V {t) = arccos = arctan ^ = - arctanr (75) 

kB k z 

between the external magnetic field and the wave- vector. 
In such a way we find the geometrical meaning of the 
dimensionless time. 

The effective friction coefficient j s (t) in the oscillator 
equation Eq. (I73p is presented by the dimensionless func- 
tion a(r) 



ls(t) = A 1s {t), 
7 s (t) ee -a(r) = 



2r 



sin 2(p. 



(76) 
(77) 



Without any restrictions we have considered k x = and 
K x = because it is related only to the choice of the 
initial time Eq. ([ilj). 

Our next step is to determine the conditions that en- 
sure small dissipation rate. Ohmic resistivity and vis- 
cosity produce an additional friction coefficient in the 
effective oscillator equation Eq. (|73")1 ; see Ref. sec. 69, 
problem 

lv = i v k + f m )k 2 (i) = Ay v , 



7„(t) = {l + T 2 )a = 



cos 2 tp ' 



v'K\ = K + v m )k 2 JA. 



(78) 
(79) 

(80) 



The upper formulae are meticulously obtained in Ap- 
pendix [DJ Taking into account the small dissipation in 
the oscillator equation (|T3")) we have to substitute the 
shear attenuation with the total attenuation 



7 = J a (f) + lv= A l = Mis + lv)- 



(81) 



The dimensionless attenuations 7 S and 7„ are more con- 
venient for further analysis of the kinetics of Alfven 
waves. 

From Eq. (]T3[) we derive an equation for the atten- 
uation of the averaged effective energy for a fictitious 
particle with a coordinate u ;- 



of the oscillator [3J, Eq. (25.5), ibid. [3j, sec. 51, problem 2 

d t E cS a -*y(<p)E eS . (83) 

Here the time dependence is implicitly included by the 
dependence of the attenuation on the angle <p, Eq. ([75]) . 

To derive the approximate oscillator equation Eq. ([72")) 
we assumed that the alteration of the wave vector com- 
ponents K(t) and respectively the variation of oc(t) (or 
7« (t) in Eq. fT3"|)) is negligible compared to the change of 
the velocity v and its derivatives and therefore all terms 
containing K x (t) and d(r) have been omitted. Such ap- 
proximation corresponds to the WKB approximation in 
quantum mechanics where we have fast-oscillating ampli- 
tude and slowly changing coefficients. The approximate 
equation Eq. (|73p will be the starting point for our fur- 
ther analysis. 

We will illustrate the WKB method for calculation of 
the attenuation on the simplest possible example of the 
system Eq. ([TO]) and Eq. (J7T]) or equation Eq. (72]). Sub- 
stitution in this equation of amplitudes oc exp(Ar) gives 
the characteristic equation 



P(A) = A 2 + K 2 



D(X) = P(X) + V(X) = 0, (84) 
d A P(A) = 2A, V(X)= ls X. (85) 



Supposing 7 S to be a small perturbative correction we 
have in the first Newtonian iteration 



A a Aq 



P(A ) = 0, 
d A P(A ) 



A = -iK z , 
1 



(86) 
(87) 
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In such a way for the energy of the wave we obtain 

E(t) oc I exp(Ar)| 2 = exp(-7 s r) = exp (-%*). 

This procedure applied to the more general system of 
equations Eqs. ([62"1 l64|) leads to the secular equation 



ls + V ' k K 2 { T )+\ -K z 

K z v> m K 2 (T)+\ 



0. 



Thus the general expression for the eigenvalues is 
A 2 + 7 A + K 2 + [ 1s K 2 (t) + v' k K\r)} v' m = 0. 



(89) 



(90) 



For high enough temperatures T 3> Ti- m the last term 
is negligible (see Eqs. (fT4l [T5l and HT))). Then we may 
ignore the slight variation of the real part of the Alfven 
frequency and account only for the imaginary correction 



A w -iK„ - -7, 



(91) 



which after a time differentiation gives the kinetic equa- 
tion for the averaged energy of the MHD wave Eq. ([83]) . 

The damping of the waves can be calculated directly 
as a ratio of the dissipated power divided by the energy. 
Initially, we have to take the real part of the wave vari- 
ables and then we have to average over the period of 
oscillations 



7u = 



(j(t,r)-E(^r)) + §E-^i((^ + ^r 
£< u 2(t,r)> + 4<B'2(t,r)> 



where the numerator represents the volume density of the 
total dissipated power Qtot and 

d 

d k = - — , (xi, x 2 , x 3 ) = (x, y, z). (93) 



Eq. ([52) is an alternative way to derive Eq. (|78p for the 
case of small shear as it is done in Appendix [Dl 



VII. KINETICS OF ALFVEN AMPLITUDES IN 
WKB APPROXIMATION 

In the spirit of quantum mechanics the number of par- 
ticles is proportional to the square of the amplitude. In 
this sense using the dimensionless amplitude of the ve- 
locity v one can introduce the mean number of "alfvons" 



n(r) = ((Rev) 2 ), 



(94) 



where the time r average is taken over the dimensionless 
period of the Alfven waves 2ttA/lua — 2n/\K z \. Accord- 
ing to the definition Eqs. (|4"6l H7|) the variables v and b 
are real, so the sign for a real part can be omitted. 

The volume density of the wave energy is propor- 
tional to the square of the amplitude and the number 
of "alfvons" 

E = 2 ^((Reti,) 2 + (Reu z ) 2 ) = pV%{\ + r 2 )n. (95) 



Here we have taken into account that the averaged ki- 
netic energy of the fluid is equal to the averaged energy 
of the magnetic field which plays the role of the "poten- 
tial" elastic energy of these transversal waves. We have 
also used the explicit form of the wave-vectors Eqs. (|55l 
l66|) . The influence of the y mode will be assessed later. In 
such an interpretation the equation for the time deriva- 
tive of the wave energy Eq. (|8"3"|) can be considered as a 
Boltzmann kinetic equation for the number of "alfvons" 



(96) 
(97) 



d T n(r) 


= -l(T)n( 


7 


= ls+lv, 




2r 




1+T 2 ' 




= + 



(98) 
(99) 



In order not to interrupt the explanation the derivation 
of the turbulence-induced source term w will be given in 
a separate subsection I VIII Bl see Eq. (|130[) . 

One can easily check that the solution of the kinetic 
equation Eq. (fTO)) reads 



n(r) 



w(t') exp 



7 (r")dr" 



dr'. (100) 



Let us first analyze in Eq. (|96[) the dissipationless 
regime of 7„ = and zero turbulence power w — 0. Using 
the integral 



exp 



ls(r"W 



cos 2 ip 



1 



(101) 



for the solution of the homogeneous linear equation we 
obtain 



n(r) ps n cos 2 ip = 



"o 



1 



(102) 



where the integration constant rig determines the spectral 
density of the waves with wave vector k parallel to the 
constant external magnetic field Bo. 

Further, in the general solution Eq. (|100|) we have to 
substitute the integral 



exp 



7,(r")dr" 



COS 2 if 



1 



COS 2 tp' 1 + T 2 

The contribution of the dissipation is given by 

1 c 



lu(r"W = 
or for the interval pointed above Eq. (|100p 



-(1 



T 2 )U- 



TU 2 U 3 



(103) 



(104) 



(105) 



where for simplicity we have introduced a shifted time 
variable 



r - r' > 0. 



(106) 
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Finally, the solution of the Boltzmann equation Eq. (|100[) 
takes the form 



1 



1 + T 2 



n(r) 



x exp < —a 



10 (t -u)(1 + (t- u) 2 ) 



(1 + T 2 )U - TU 2 + -U 3 



du. (107) 



In order to evaluate this integral we will assume wave- 
vector independence of the random turbulent noise 
w(K(t)) w const. For the real physics of accretion disks 
the bare viscosity is evanescent and we have to analyze 
the above integral on the u variable for very small values 
of a <C 1. This means that for values of |r| of the order 
of 1 we have to take into account very large values of the 
u 3> 1 variable. For large u in the integrant we have to 
make the approximations 



1 + {t-u) 2 ra u 2 , 



(1 + T 2 )U - TU 2 + -U 3 

O 



1 3 

-u J 
3 



(108) 
(109) 



and for the integral in the above Eq. (|107|) we derive 
according to Eq. ()102[) the r-independent evaluation 



w a 



(110) 



In such a way for the low dissipation limit we derive an 
explicit expression for the number of "alfvons" propagat- 
ing along the magnetic field 



Aw 



n 



{vk + v m )k 2 



(111) 



The last term is expressed via the physical variables of 
the initial problem and this is the central result of our 
analysis of the kinetics of slow magnetosonic (Alfven) 
waves. Later, we will perform statistical averaging of the 
dissipated power, but, before that, in the next section, 
we will analyze in short the influence of weak turbulence 
on the MHD equations. 



VIII. WHITE NOISE TURBULENCE 

A. Langevin-Burgers noise for Alfven waves 

In Langevin-Burgers approach to the turbulence the 
random noise is described as a white one for which the 
following correlation may be used 

(F(ti, n)F(t 2 , r 2 )) = Tp 2 6(h - t 2 )£(n - r 2 )l, (112) 

The dimension of T shows that this correlation is conve- 
nient not for forces, but rather for accelerations 



f = (acceleration) x (volume) x (time). (113) 



We make a Fourier transform of the acceleration 

F(i,r) = £)F fc (t)exp(ik.r), (114) 

k 

F fc (t) = y"F(t,r)ex P Hk-r)^ (115) 

for which we apply periodic boundary conditions 

F(t,x + L)=F(t,x) (116) 

where L is the characteristic length of the system V = L 3 . 
This sets the following restrictions for the wave vectors 

(k x , ky, k z ,) — —j-(n x , n,y, n z ), (H7) 

where the numbers can take only integer values 

n x , %, n z = 0, ±1, ±2, ±3, ... (118) 

In this case the difference between two wave vectors is 
Afc = 2tt/L and for great lengths the sum turns into an 
integral 



1 E 



d 3 fc 

(2^)3 



(119) 



The correlation for the coefficients in the Fourier series 
is also a white noise 



{F;(h)F k (t 2 )) = t^6(h - t 2 )S pk 



(120) 



where S p k is the symbol of Kroneckcr. Taking into ac- 
count the upper relation for the x component of the ac- 
celeration we have 

(Fl x (h)F k , x (t 2 )) = '^S(t 1 -t 2 ). (121) 

In our work we are interested in the presence of a white 
noise in the MHD system (|62|) - ([64]) causing the pri- 
mary birth of the Alfven waves. That is why we want to 
know the correlation for the dimensionless density of the 
external force 



where 



(r(n)f(T 2 )) = T6(n-r 2 ), 



/WeM r r 



pAV A 



VlAV 



(122) 



(123) 



and the parameter T is dimensionless. 

This dimensionless random force has to be taken into 
account as an external force in the right hand side of 
Eq. ([70]) which now takes the form 

d T v x = a(r) v x - K z b x + /(t) (124) 

and analogously Eq. (fT2"|) reads 

v = a(r)v - K 2 v + /(r). (125) 
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Introducing the dimensionless displacement 

x{t) = ( v{t')At' (126) 

J — oo 

we obtain 

x = - 1s {t)x-K 2 z x + f{r). (127) 

In such a way we arrived at the necessity to analyze 
the behavior of classical harmonic oscillator with time- 
dependent friction and a white noise as an external ran- 
dom force. 



B. General theorem for stochastic heating and 
income term in the effective Boltzmann equation 

The special cases considered in Appendix|E]for a Brow- 
nian particle Eq. (|E15|) . an oscillator without friction 
Eq. (jE3T|) and a free particle Eq. (jE36]) give us a hint 
that there is a general formula for a white-noise stochas- 
tic heating. Let us apply this common result to the effec- 
tive oscillator equation Eq. (|127|) for the displacement of 
the plasma x. In case of negligible 7 S Eq. (|E31[) implies 



I 

2' 



(128) 



where the dimensionless Burgers parameter T is defined 
in Eq. (|123[) . At slow heating the virial theorem (|i 2 ) = 
(^K 2 x 2 ) gives that 



d T <± 2 ) = ^-. 



(129) 



According to the determination Eq. (|126[) x — v and 
Eq. ([94| reads n(r) = (± 2 ) . In this way, in case of zero 
friction 7 = 0, we obtain 



d r n(r) 



r 

2 



2AVIV 



(130) 



which is the source term in the right hand side of the 
Boltzmann equation Eq. (I96| . The turbulent random 
forces generate magnetosonic waves. Figuratively we 
may say that "alfvons" are born from the turbulent sea 
foam. The other terms in the master equation Eq. (|96j) 
describe the dissipative decay with rate 7„ and the "las- 
ing" at 7 S < 0. 



IX. PERSPECTIVES AND CONCLUSIONS 

The discussed mechanism requires simultaneous pres- 
ence of both magnetic field and turbulence - ingredients 
which arc pointed out practically in all present works on 
the problem. As we have analyzed MHD waves in an 
incompressible fluid the gas pressure turned out to be ir- 
relevant to the shear stress. Therefore we consider that 



in the Shakura-Sunyaev phenomenology the gas pressure 
p should be replaced with the magnetic one p B . Thus for 
the momentum transport we have 

<Jtpr — a^p B for disk, a zr — a u p B for column, (131) 

which is a rather small correction to the Shakura-Sunyaev 
hypotheses . 19 ' 20 The detailed calculation of the dimen- 
sionless parameter a u requires complete analysis of the 
kinetics for all modes in the shear flow of the magnetized 
plasma. This is undoubtedly a very complicated task and 
it is worth making a qualitative assessment for the origin 
of the effective viscosity in a plasma with evanescent ini- 
tial one. The most important characteristic of the shear 
flow in a magnetized plasma is the "lasing" phase when 
the wave draws energy from the shear flow. Without 
dissipation this increment reaches gigantic scales^ 

Complete analysis of the system, of course, requires de- 
tailed numerical integration. Despite of this we can give 
a qualitative explanation that the plasma heating and 
the arising of a huge effective viscosity is caused by the 
"laser" amplification of the convective instability. This is 
a rather universal mechanism playing a significant role in 
our solar system in the diminution of the angular momen- 
tum of our sun, but also it is the mechanism responsible 
for the gigantic energy produced in quasars and AGNs. 
In this way we offer a qualitative answer to the question 
why do the most powerful sources of light in the universe 
shine - namely, because of the instability of the slow mag- 
netosonic waves in the shear flow. We wish to point out 
that our theory for heating of accretion disks is a conven- 
tional one, which does not require additional hypotheses, 
but only numerical calculations. The energy is produced 
in the bulk of the disk but finally emitted to its sides. 

The analyzed geometry B = (0, 0, B ) corresponds, 
for example, to an accretion column above the magnetic 
poles of a neutron star. However, we suppose that the 
convective instability of the slow magnetosonic waves is a 
general property of the magnetized plasma with a shear 
flow. We believe that an arbitrary orientation between 
the shear flow and the magnetic field would lead only to a 
dimensionless multiplier of the order of unity in our final 
result. 

The used Langevin-Burgers' approach is a first approx- 
imation for a self-consistent treatment of the turbulence. 
We can modify it including, for example, a wave-vector 
dependence of the noise functions T(k) in order to simu- 
late the spectral density of the velocity pulsations. Some 
parameter related to the spectral density of the turbu- 
lence will indispensably participate in the final result for 
the stress tensor. Moreover we have to include the self- 
sustained turbulence in a realistic global model for the 
accretion disks or accretion streams. In this direction we 
are not giving a pret-a-porter prescription. The purpose 
of our work is only to give an idea how the turbulence, 
magnetic field, and shear flow play together in the most 
powerful sources of light in the universe (for the mecha- 
nism of this shining we do not know more than the girl- 
friend of H. Bete on the energy source of stars, cf. the 
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story told by Feynman in his famous lectures). Simul- 
taneously we can observe traces of a big viscosity of the 
protoplanetary disk in the angular momentum distribu- 
tion in our solar system X We just wish to rise the corner 
of the curtain and see the regally play by: 1) the tur- 
bulence 2) magnetic field and 3) shear flow creating the 
most efficient engine in the universe. 

The work of this engine is related to the kinetics of the 
spectral density of slow magnetosonic waves in a magne- 
tized turbulent plasma with a shear flow. The heating, 
effective viscosity and stress tensors are statistical con- 
sequences of this spectral density. We advocate that the 
"lasing" of "alfvons" is a key detail in the accretion of 
many compact astrophysical objects in order to observe 
the universe in its present form. 

The last problem which we wish to speculate on is 
the applicability of the Burgers' approach to magneto- 
hydrodynamics. According to the Lighthill theory^ the 
intensity of the emitted sound is proportional to the ra- 
tio of the velocity pulsation and the sound velocity to a 
high power (v/v s ) 5 , see also Ref. Q, sec. 75. Qualitative 
considerations give that the velocity of sound v s has to 
be substituted by the Alfven velocity Va -C v s for small 
magnetic fields p B <C p. This leads to the conclusion 
that the transformation of the turbulent energy to Alfven 
waves can be very effective and those energies could be 
comparable. This qualitative property was pointed out 
for the physics of solar plasma long time ago, see for ex- 



ample Refs. [2lll22l ]. In such a way in order to model the 
influence of the turbulence on MHD waves we have to 
use high values of the noise intensity T. This parameter 
can be determined in such a way that the spectral den- 
sity of the waves' energy becomes equal to the spectral 
density of the turbulent energy at a given cut-off wave- 
vector k c . For qualitative computer simulations in the 
inertial regime we can use these typical cut-off values as 
initial conditions and solve the MHD equations without 
the stochastic force. This could be useful for a numerical 
calculation of the spectral density of Alfven waves, i.e. 
"momentum" distribution of alfvons in our conditional 
terminology. 

We conclude that it is necessary to revise the theoreti- 
cal models of disk accretion and we might expect appear- 
ance of a new direction in the theoretical astrophysics, 
incorporating the methods of statistical physics as a key 
detail in global magnetohydrodynamic models for forma- 
tion of compact astrophysical objects. 
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APPENDIX A: MATRIX PRESENTATION FOR 
LAGRANGE-EULER TRANSFORMATIONS 

The transition between Lagrange and Eulerian coor- 
dinates for the space- and wave-vectors is essential for 
the present paper. For convenience of the reader we are 



giving these relations in a transparent self-explainable 
matrix form: 

(k|r) = (k|r), (k| = (k k y ~k z ), (Al) 



\k x ky k z j 



1 CP 
(k x k y k z ) | 1 
,r0 1, 





|k(t)> = 




(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



APPENDIX B: DYNAMICS EQUATIONS 

Let us consider the movement of incompressible mag- 
netized homogeneous plasma. The time evolution is given 
by the system^ 




= 1, 










f 1 







• 


1 


T i 





1 











Vo 





r 



p(<9 t v + v- Vv) 



-Vp- — (B x rotB) 
Mo 



9(B+(vVB) = (B-V)v+ AB. 



?7Av + F, 
(Bl) 

(B2) 



We presume small variations from the stable laminar cur- 
rent (v = U +u , B = B +B , p = po+p )■ Then, using 
the assumptions for the direction of the external mag- 
netic field and the shear flow velocity (Bo = Boe z , Uo = 
Axe z ), we find 

(v-V)v=(U +u)-V(U + u) 
= (U • V)U + (u • V)U + (U • V)u 

= Au x e z + Axd z u. (B3) 

Analogously, for the mixed multipliers we obtain 

(v • V)B = (U + u) • V(B + B') 
= (Uo • V)B + (u • V)B + (U • V)B' 

= Axd z B, (B4) 
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and 



we obtain a system of ordinary differential equations 



(B ■ V)v = (B + B ) ■ V(U + u) 
= (Bo ■ V)U + (B' ■ V)U + (B ■ V)u 

= B d z u + AB x e z , (B5) 
Vp = Vp' , AB = AB' , Av = Au. (B6) 

The time derivative of both variables is expressed by their 
varying components 



<9 t v = <9 f u, d t B = d t B . 



(B7) 



Having in mind the above mentioned assumptions the 
vector product of the magnetic field and its rotation may 
be rewritten in the form 

B x rot B = B x rot B' 

'(B ) a (rotB% - (Bo) z (rotB'y 
(B ) 2 (rotB% - (Bo)x(rotB' 
.(BoMrotB'^ - (BoyrotB')^ 

r -B (rotB') y \ (d x B z - d z B x s 



= \ ,(rotB ;, J = B \d y B z -d z B y 

Consequently, Eqs. (|B1[ IB2[) turn into 

p(d t u + Axd z u + Au x e z ) = — Vp 



(B8) 



So 



d z B — d x B z 



1 



(B9) 



d t B + Axd z B = B d z u + AB x e z + — AB . (BIO) 

per 



APPENDIX C: DERIVATION OF THE MHD SET 
OF EQUATIONS WITH RANDOM NOISE 

Substituting a plane anzatz waves in the linearized 
MHD Eqs. HUD 

u(t,r) = -iV A v(T)exp(ik(t)-r), (CI) 

B'(t,r) = S b(r) exp(ik(t) -r), (C2) 

p(t, r) = pVl P(t) exp(tk(f) • r), (C3) 

F(t,r) = -ipAV A {'{T)exp(ik{t) -r), (C4) 

where according to Eqs. (f3"o) . and (|4*Tj) 

K • b = 0, K • v = 0, (C5) 



k(t) = yrK(r), K(r) 




K 



K z e z , 



(C6) 



v = -v x e z + KP 

K z b x — K x b z 

K z b y -K y b z ) - i4K 2 r f. (C7) 



b = K z v + b x e z - v' m K 2 h. 



(C8) 



Here r = At, the dot operation stands for r- 
differentiation d r and we have introduced dimensionless 
notations for both kinematic and magnetic viscosities 



A 



A 



(C9) 



All dimensionless variables in the upper system are ex- 
pressed in specific units: for velocity Va, pressure pV\i 
magnetic field Bq, time 1/A, wave- vector A/Va, kine- 
matic viscosity V%/A, acceleration AVa, density of force 
pAV A , and length V A /A. 

Time differentiation of Eq. (|C5j) gives 

K x v x + KyV y + K z v z - K z v x = 0. (CIO) 

The substitution here of v from Eq. (|C7[) gives for the 
pressure 

P=-bz + 2^v x - j^(K x f x + K y f y + K z f z ), (Cll) 

and back substitution of the dimensionless pressure in the 
x- and y-components of Eq. (|C7|) gives the final equations 
for the acceleration 

K K 

= 2 * K 2 X Vx ~ Kzbx ~ v K Vx 

+f' x - ^(K x f x + K y f y + K z f z ), (C12) 
K K 

iL = 2- * " »■ ••'^ 2 - 



K 2 



) x - K z b y - V K Vy 



K, 



+fy ~ j£&*fL + Kyf y + K z f z ). (C13) 

For the magnetic field we take the x- and y-components 
from Eq. 



b x = K z v x - v m K b x , 
by = K z v y -v' m K 2 b y . 



(C14) 
(C15) 



The corresponding equations for the z-componcnts 
v z = (-l + 2^v x -K z b z -v' k K 2 v z 

+fz - ff + K vfy + K*f' z ), (C16) 

b z = K z v z +b x -v' m K 2 b z (C17) 

are not necessary to be solved; the solutions are given by 
Eqs. jC5} 

v z = -(K x v x + K y v y )/K z , (C18) 
b z = -(K x b x + K y by)/K z . (C19) 
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1. 2D restriction with K y = 

In the special case of K y = we have separation of 
variables and the y-component of the acceleration is 



Vy = -K Z by - V K Vy + f' y 



(C20) 



where for the Langevin force we suppose a white noise 
correlator 



(f'(n)f / (T 2 )) = ir<5(Ti-T 2 ), r 



V 2 AV 



(C21) 



In such a way we derive Eqs. 

In the present work we will use the K y = case for 
a model evaluation of the statistical properties of Alfven 
waves. 



Test example of short wavelengths and small 
attenuation 



As a test example let us finally analyze the textbook's 
case of a negligible viscosity and big enough wave-vector 
K z . The systems for v x , b x and v y , b y have approximate 
time dependent solution oc exp (— z|if^|r) . This means 
that both x- and y-polarized modes have dispersion of 
Alfven waves u> ~ A\K Z \ or more precisely 

w(k) =uj a = V A \k z \ - % -vk 2 , uj A -»A : vk 2 . (C22) 

This imaginary term may be obtained as the first Newton 
correction in the method described above (see Eq. (|87|0 . 
Let us mention that the attenuation given by the small 
imaginary part of the frequency is an isotropic function 
of the wave vector, cf. Ref. 0, sec. 69, problem. 

We analyze an inviscid approximation and that is why 
for both MHD branches the oscillations of the velocity u 
are transversal to the wave-vector k. For the y-mode the 
velocity oscillation u is perpendicular to the k-Bo-plane; 
this is the true Alfven wave. For the x-mode, which we 
are analyzing in the current paper, the displacement of 
the fluid and the magnetic force lines lies in the k-Bo- 
plane. For finite compressibility the x-mode is hybridized 
with the sound, that is why it is often called slow mag- 
netosonic wave even if at low magnetic fields p B <C p 
its dispersion is given by the Alfven waves' dispersion 
u A nV A \k-B \/B . 



APPENDIX D: ATTENUATION COEFFICIENT 
FOR ALFVEN WAVES 

In order to derive the dissipative function for Alfven 
waves spreading in a highly conducting plasma, we have 
to know the energy and its variation in time due to kine- 
matic and magnetic friction. In this section we are in- 
terested only in the average dissipation, that is why it 



is convenient to work with averaged (with respect to the 
wave's phase) quantities. The density of electromagnetic 
energy is given by 



(Dl) 



where we have used that for a highly conducting media 
the electromagnetic energy E 2 /2sq cx 1/ct 2 is negligible 
compared to the magnetic one. 

We are interested in the wave part of both the mag- 
netic and the velocity fields, as the constant part does 
not play any role in the dynamics of the system. Rewrit- 
ten in terms of the dimensionless wave components (see 
Eq. (|47|) ) the expression for the energy density looks like 



b 2 p E 



(D2) 



where the magnetic pressure p B — Bq/2u,q has already 
been introduced in the text above (sec. HIIi Eq. (JTJ) ) . Ac- 
cording to the equipartition (Alfven) theorem the Alfven 
waves kinetic energy density equals the magnetic one 



£kin = -^pV 2 = v 2 p B . 
Then having in mind that 



(D3) 



(D4) 



for the averaged volume density of the total energy we 
obtain 



2P B [(1 + r 2 ) (vl) + (v 2 y )] 



(D5) 



For an estimation of the attenuation coefficient one also 
needs to know the averaged with respect to the waves' 
period volume density of the total power Q tot , which is 
a sum of the averaged kinetic and magnetic one 



Qt 



ikin 



(D6) 



The kinetic power is simply the time derivative of the 
kinetic energy and since we are interested in the wave 
component of the velocity it is given by 



(D7) 



where the separation of variables Eq. (f3"Tj) and the lin- 
earization Eqs. (|3"T1 13"9"|) have been used. 

Rewritten in the dimensionless variables Eq. with 
the incompressibility condition Eq. (|52[) applied this 
reads 



+ Vy 1 + 



Qkin = -2Ap B v ■ d T v 

K K 

° K 2 Z '" X ~ Kzbx ~ V 'k K2v x + 9 fx 



K K 

V K2 Z v x - K z b y - v' k K 2 v y + hfy 



(D8) 
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When we take an average value of the upper expression 
only the quadratic terms remain 



(by 



0. 



(D9) 



J yi — \ u vi 

therefore for the kinetic part of the average power we 
have 



(Qkin) = ~2Ap E 



1 



K r K z 




(D10) 

If we use the relation Eq. (|55ll for the special case when 
K y = this yields 

(Q kin > = 2Ap B K 2 ((1 + r 2 )\' k {v 2 x ) + (1 + r» 2 )) 

+2A Pb t(v 2 x ). 

(Dll) 

Now we have to calculate the Ohmic part of the power 
in the magnetostatic approximation Eq. (jTTJ) 

,2 



Qc 



(rotB) z 



2p B ^ m (rotb) 2 = 2A P y m {K x b) 2 



(D12) 



where E' is the effective electric field in the Ohm's law 
Eq. (TTU|) . b is the dimensionless magnetic field defined in 
Eq. (143)) and v' m is the dimensionless magnetic viscosity 
given in Eq. (T57|) . 

In the particular case when K y = the upper expres- 
sion in terms of the dimensionless time r turns into 

Q cm - 2Ap B v> m K 2 (l + r 2 ) [(1 + r 2 )b 2 x + &£] . (D13) 

Therefore if we take into account Eqs. (JdH EH |DTT|) 
and the equation above for the average density of the 
total power dissipated in the fluid we have 

Q tot = 2A Vb v'K 2 X\ + r 2 )((l + t 2 )v 2 x + vl) + 2A Pb tv 2 x , 

(D14) 

where v is the total viscosity defined in Eq. ({551) . 

Now it can be easily shown that the attenuation coef- 
ficient takes the following form 

A(l + T 2 yK 2 + -^. (D15) 



_ Qtot 

7 : 



1 



In this way we derived the approximative (in case of small 
dissipation) equation Eq. ([81]) . 



even have not heard about Langevin-Burgers' approach 
to turbulence^ 3 .. That is why instead of referring to text- 
books far away from our current problem 26 we will give a 
pedestrian introduction of all the necessary basic notions. 
In the framework of the used notations we will introduce 
the necessary mathematics. Our first illustration will be 
the diffusion of a Brownian particle. The theory is anal- 
ogous to the Nyquist^i theory for the thermal noise in 
electric circuites. 

The literature on "Burgers turbulence" is really 
huge http://google.com search gives 29, 000 items 
for 2007, or if we introduce a google index func- 
tion it will read: Q ("Burgers turbulence", 2007) = 
29, 000. Analogously, for Langcvin turbulence we have 
g (Langevin turbulence, 2007) = 113, 000. We will men- 
tion only several reviews on this subject^. 

Our goal in this subsection is to derive an explicit and 
physically grounded expression for the random noise in 
the Boltzmann equation Eq. (|M|) . We shall start our 
deduction by consideration of the kinetic equation for 
a Brownian particle, moving with velocity v under the 
influence of a random force f (t) across a given medium 
with a "friction coefficient" 7 



mdtv = —7717V + f (i). 



(El) 



For the sake of simplicity from now on we shall examine 
the propagation of the particle in x direction, so that we 
shall be interested only in the projections of the velocity 
and the external force along this axis, which we shall 
designate with v and /. The upper ordinary differential 
equation may be easily solved with the help of the Euler 
method 



v(t) = C(t)exp(-jt). 



(E2) 



We substitute the velocity from Eq. (|E1[) with the given 
expression and obtain another ordinary differential equa- 
tion with separable variables for the unknown function 
C(t) 



exp(-7i)d t C(t) 



fit) 



whose solution is 



C(t)=vo+ f ^exp(7ti)dti. 

Jto 



III 



(E3) 



(E4) 



APPENDIX E: STOCHASTIC HEATING 



The constant of integration vq here stands for the initial 
velocity. The explicit form of the velocity becomes 



1. Kinetic equation for the kinetic energy of a 
Brownian particle 

The Langevin 2 ^ approach for a treatment of stochas- 
tic differential equations is practically unknown in astro- 
physics; most of the actively working in this field people 



v(t) = exp(-Tt) I vo + £ exp( ~ / , !< 1/ , j . ( Er, ) 



■/?) 



Now let us average the square of this velocity. For this 
purpose we need to know the correlation of the noise. Fol- 
lowing the classical works by Langevin and Burger o 28 i 33 
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we suppose that the external force correlator has the sim- 
plest possible form of a white noise 



f(h)f(t 2 )) =T6(h-t2). 



(E6) 



We take into consideration that the random force has 
zero mean value (f(t)\ = and for the average of the 
square of the velocity we obtain 

(v 2 (t)) = exp(-2 7 i) 
x K+ / / rS ^ h ~ ^ exp(7(ti + f a ))dfidf a ) 



to J to 



m- 



exp(-2 7 i) [ v 2 + f i exp(27ii)d<i 
V Jto m 



(E7) 



This noise averaging is the most important ingredient in 
the present derivation. It reduces the mechanical prob- 
lems to a statistical problem solvable by the Boltzmann 
equation. In the next subsection we apply this approach 
to a harmonic oscillator and later on to the amplitude of 
magnetosonic waves. 

We are now ready to calculate the mean kinetic energy 
of the Brownian particle 

(E kin (t)) = ^(v 2 (t)) 

mv 2 f /"* 

= — ^ exp(-2 7 i) + — exp(-2 7 t) / exp(2 7 t 1 )dfi 
2 2m J to 

= ^0 exp(-2 7 t) + {l _ exp(-2 7 (i - to))} (E8) 
2 4m 7 

To rewrite this expression in a more convenient way we 
introduce the initial kinetic energy Eq — Ekmtyo) & n( i the 
equilibrium one E = T/4mj. Thus the average kinetic 
energy reads 

(E kin (t)) = E exp (-l^j+E^l- cxp (~^^ 

E . (E9) 

where te = l/2 7 corresponds to the relaxation time in 
the atomic physics and T is the analogue to the width 
of the Lorentz function. Time differentiation of Eq. (|E9|) 
gives the kinetic equation for the kinetic energy 



d t (E kin (t)) = -— [E cxp [- — 



Sexp 



(t - t Q ) 



(E10) 



In this way we came to the well-known Boltzmann equa- 
tion for the variation of the average kinetic energy 

£ (E kin (t)) = -— {(E kin (t)) - E) . (Ell) 
d* te 



The first term here is responsible for the energy expendi- 
ture and consequently is set by the dissipation function, 
whereas the second stands for the energy income (i.e. a 
positive power), caused by the effect of the random noise. 

Now let us make a short analysis of the result. In 
the case of equilibrium (after several relaxation times, 
t — > oo) the mean value of the kinetic energy equals half 
the absolute temperature which we will designate with & 



=le=E=-L. 

2 4m 7 



(E12) 
(E13) 



Hereby, we can deduce the connection between the dissi- 
pation coefficient 7 and the fluctuations T 



T = 2m8 7 , 



(E14) 



which is a special case of the fluctuation-dissipation the- 
orem. 

Our idea is to apply notions from the equilibrium 
statistics in the non-equilibrium case, therefore we shall 
use 8 as a white noise parameter in stead of T (a non- 
equilibrium analogue to the equilibrium temperature). In 
terms of the averaged power of energy dissipation Q, the 
fluctuation-dissipation theorem reads 



Te 



7 e 



2 m 



(E15) 



The first special case of the fluctuation dissipation the- 
orem is the relation between diffusion coefficient and the 
mobility of a Brownian particle. We will remind some 
details. Let us consider the motion of a Brownian parti- 
cle as a diffusion, for which the role of the concentration 
is played by the probability to find the particle at a given 
place in the volume of the fluid. The diffusion equation 
is 



dtn(x, t) — DAn(x, t), 



(E16) 



where n(x, t) stands for the concentration. The square 
of the average distance, which the particle passes for an 
interval of time t is given by Ref. Q, sec. (60) 



(E17) 



According to the Sutherland-Einstein relatio n 27 ' 29 i 30 the 
diffusion coefficient D is proportional to the mobility 
/x and the temperature in the equilibrium statistical 
physics, which under our assumption may be expressed 
as 



D = fj,Q. 



(E18) 



The mobility may be derived from the equation of motion 
of the particle, roaming under the influence of a time 
independent random force f 



mdtv 



-m 7 v - 



(E19) 
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In the stationary case the driving force has to balance 
the friction 



Thus for the correlator we have 



Vdr = M/- 
Hereby for the mobility we obtain 

1 



and the diffusion coefficient takes the form 

e 



D 



(E20) 
(E21) 

(E22) 



Our next step is to consider the equation of motion for a 
harmonic oscillator under a random noise. 



2. Oscillator under a white noise 

In a self-consistent linearized approximation the prob- 
lem for the wave propagation is reduced to independent 
oscillator problems for all wave vectors. That is why the 
Brownian motion of a harmonic oscillator is a key detail 
of our statistical theory for the spectral density of mag- 
netosonic waves. We are starting with the equation of 
motion for a harmonic oscillator with an external white 



-miv 2 x 



•f(t). 



(E23) 



We assume classical oscillator (i.e. [x, p] =0). It is con- 
venient to examine the equation of motion in the phase 
space. For that purpose we introduce complex variables 



c = mux + ip 
c* = mux — ip, 



(E24) 



where p = mi is the particle momentum. 

The energy of the oscillator (proportional to the num- 
ber of particles) can be expressed via these variables 



E = -^-c*c 
2m 



V 

2m 



1 2 2 

—muj x 
2 



(E25) 



We are interested in the explicit form of the mean power 
Q = d t (E(t)) . That is why we want to know the correla- 
tion between c and c* . In order to find it first we have to 
determine their explicit forms as functions of frequencies, 
random noises and time. The time derivative of c looks 
like 



d,c 



mux + ip = Lup - 



f(t) 



= — iuj(ip + mu>x) = —iu>c{t) + if it). (E26) 



For this ordinary differential equation we seek a solution 
in the form of harmonic oscillations 

c(t) = I Cq +i / f(ti) exp(iwti)dii I exp(— iuit) 



c*(t) = [C -i J ffa) exp(-iwi 2 )dt 2 J exp(— itot). 

(E27) 



(c*(t)c(t)) = \C Q \ 2 
/(*i)/(*a) ) cx P M<! - h))dt x &b 2 . (E28) 



We consider a white noise for which (fitijfiti) 

T5(ti - i 2 ). Then Eq. (|E28|) turns into 

(c*(t)c(t)) = \C \ 2 + T(t-t Q ). (E29) 

This averaging is one of the most important details of 
the present theory - it reduces the wave problem to a 
statistical one solvable by the Boltzmann equation. 

Finally for the averaged energy of the oscillator under 
the random noise we have 



(E(t))=E + —(t-t ), 



(E30) 



where the initial energy is E = ^ \C Q \ 2 . The produced 
power is the same as in Eq. (|E15|) 



Q = — . 

2m 



(E31) 



Here we wish to mention that analogous scenario of lin- 
ear dependence of the energy of ocean waves driven by 
turbulent fluctuations of the pressure was proposed by 
Jeffries, Fillips, Feynman and Hibbs. This model, how- 
ever, could be applicable only before the creation of a 
system of parallel vortices which are the real interme- 
diary between the wind and waves (l2lfl3j : see also the 
references therein. 

The nature of the friction force and the external poten- 
tial is irrelevant to this result because by definition the 
white noise is a very fast process. In order to illustrate 
this general theorem we will analyze in the next section 
the averaged power of a free particle under a white noise. 



3. Stochastic heating of a free particle 

This phenomenon is analogous to the Fermi model for 
acceleration of cosmic particles. The equation of motion 
for a free particle reads 



= /(*), 



(E32) 



where fit) stands for the random (under our considera- 
tion white) noise. Hereby, after integration, for the ve- 
locity vector we find 

"CO = / — dt i + v o, (E33) 
J to m 

where vq is the initial velocity. Using this explicit expres- 
sion the average kinetic energy reads as 

m)) = ^{At)) 
= 'LS* f t ( h{t) ^ {2) ) dtidt2 + y ^ ■ (E34) 
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Here we apply the correlation of the white noise 
A(*)/2(2)y = tS(ti — t 2 ) and integrate with respect 
to time. Then the average kinetic energy becomes 

(E(t))=Eo + 1 ). (E35) 

2m 

Finally a time differentiation gives the power of the white 
noise acting on the free particle 

Q = d t (E(t)) = ^-. (E36) 
2m 

Again this is the same result as in Eqs. (|E15|) and (|E31[) . 
This constant power is applicable if the typical relaxation 
times of the particle are much bigger than the character- 
istic times of the noise and we can approximate it as 
a white one. When we apply this result to the averaged 
amplitude of the magnetosonic waves this constant power 
gives the constant income term w in the Boltzmann equa- 
tion for the spectral density of the MHD waves. 

APPENDIX F: NUMERICAL ANALYSIS 

1. Pure hydrodynamics 

Our first step in the numerical analysis is to investi- 
gate the influence of the shear flow for amplification of 
initial perturbation in framework of pure hydrodynamics 
at zero magnetic field b = 0. For an ideal fluid v' k = 
the equation Eq. for K y = case takes the form 

d T v x = a{r)v x (Fl) 

and have the solution 

v x (r) = — l —, (F2) 

1 + T Z 

which is depicted at Fig. 1. This time dependence is 
common for all K z for which the dissipation is negligible. 



2. Nonzero magnetic field. Analytical solutions 

We will perform analysis of the model case for K y = 
0. The y-system Eq. (|63|) and Eq. ([65)1 describes Alfven 
waves for which the shear have negligible influence. The 
shear flow is important for the slow magnetosonic waves 
for which the velocity oscillations and variations of the 
magnetic field are in the plane of the wave-vector and 
the constant external magnetic filed. Only for this waves 
we have significant amplification by the shear flow which 
we will investigate in the beginning for an ideal fluid for 
which the x-system Eq. ([62]) and Eq. ([64]) reads 

d T v x = a(r)v x - K z b x (F3) 
d T b x = K z v x . (F4) 



> 
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FIG. 1: Time dependence of velocity in a shear flow at zero 
magnetic field; pure hydrodynamic solution. 
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FIG. 2: Energy of an Alfven wave as function of dimension- 
less time r. An Alfven wave even at slow viscous damping 
attenuates heating the plasma and new one is generated by 
the turbulence. We the born after billions years see the light 
of quasars created by amplification of MHD waves by shear 
flow. Parameters of the example: K z = 0.7, v' m = 0.00001, 
v' k = 0.00001, r € (-30, 120), 6(-30) = 0, and t>(-30) = 1. 



5000 - 



4000 - 



H 3000 - 



2000 - 



1000 - 




1 1 1 1 1 1 1 1 1 

-40 -20 20 40 60 80 100 120 

T 

FIG. 3: "Lasing of alfvons" as at zero damping. Ampli- 
fication of energy of the slow magnetosonic wave by shear 
flow. Parameters: the same as Fig. 2 but without friction 

Vm = V k = 0. 
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FIG. 4: "Lasing of alfvons": the time dependence of the ficti- 
tious Schrodinger equation amplitude t[>(r) oc B x {r)/\/l + t 2 
for: K z = Vol, ^(-150) = 1, d T ^(-150) = 0. 




FIG. 5: Effective potential 1/(1 + t 2 ) 2 and effective energy 
Kz —0.1 for the solution of Schrodinger equation depicted at 
Fig. 4 . 



The substitution of d T v x from Eq. (1F3|) in time differen 
tiated of Eq. jF4} gives 

d 2 T b x + a(T)d T b x + K% = 0. 



(F5) 



With the help of the substitution 



6 x (r)=V(r)Vl + r 2 



(F6) 

we arrive at the effective Schrodinger equation for the 
alfvon amplitude 

. f vl 1 \ „/. _ -2 ■ 2m 



fl+T 



(F7) 

For small enough effective energies K 2 we have a clas- 
sical forbidden region where the magnetosonic waves are 
amplified. 

This equation has two linearly independent solutions 
even (g) and odd (u) 



V g (0) = 1, d T V> g (0)=0, 

V> u (o) = o, d T ^ u (o) = i, 



(F8) 
(F9) 



FIG. 6: Even wave eigenfunction: i/> g (0) = 1, d T ip g (0) = 0. 
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FIG. 7: Odd wave eigenfunction: ^u(0) = 0, d T ^ u (0) = 1. 



^ g = VT+^H^O, -\, 0, - l±p-; -r 2 ), (F10) 



^u=r V / l+^H c (0,+i 0,- 



Kl l + Kl 



4 '' 



-r 2 ),(Fll) 



and the general solution we will represent by the linear 
combination 



(F12) 



The confluent Heun function H c (z) obeys the differen- 
tial equation 



z z — 1 



with initial conditions 



y 



Apaz — a 



V = 0, (F13) 



2/(0) = 0, y'(0) = - (F14) 

7 

and close to z — has the Taylor expansion 

oo 

E c (p,a,-y,6,a;z) = J2 C nZ n i (F15) 

n=0 
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where for the coefficients we have the initial conditions 
C_i = and Cq = 1 and recursion 

C n+1 = -(g n C n + h n C n -x)/f n , n = 0,1,2,... (F16) 

where 

g n = n(n-4p + j + 5-l)-(T, (F17) 

h n = Ap{n + a- 1), (F18) 

fn = -(n + l)(n + 7 ). (F19) 

For large enough arguments we have the expansion 

oo 

H c (p, ai7 ,^;z) = ^Ci°=' Z a - (F20) 

n=0 

where only the recursion functions are different 



7 M 



(a + n)(n-Ap-j-S+l)-a, (F21) 



= -(a + n-l)(a + n-7), 
- -4p(n+l). 

3. Simple quantum mechanical problem 



(F22) 
(F23) 



In order better to analyze the effective MHD equa- 
tion Eq. (|F7|) we will solve the corresponding quantum 
mechanical problem when we have tunneling trough the 
barrier 2mU/H 2 = 1/(1 + t 2 ) 2 , supposing that if) is a 
complex function. We have falling wave with unit ampli- 
tude, reflected wave with amplitude R and transmitted 
wave with amplitude T 



ip(r — ► — oo) 
%l)(r — > +oo) 



exp(iK z r) + Rexp{~iK z r), (F24) 



Texp{+iK zT ). 
Using the asymptotic of the eigenfunctions 

D g cos(K z t — </> g ), for t — > — oo, 



D g cos(K z t + </) g ), for r — > +oo, 



D u cos(-ftr z T 



for r 



-00, 

-oo 



(F25) 



(F26) 



(F27) 



D u cos(_ft' z r + U ), for r 
and solving the matrix problem we obtain 

C^Dg = exp^ g ), C^D U = - exp(# u ) (F28) 
Then for the tunneling coefficient we derive 
V = \T\ 2 = s 2 g , s ug = sin(0 u - g 
Let us introduct notation 

Vu = 4>u 

then Q gu = <f> g - (p u . 
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FIG. 8: Coefficient of transmission through quantum barrier 
as function of wave- vector D — cos 2 Qug- This coefficient 
parameterized amplification of alfvons A = 2/D — 1. 



4. MHD and real ip 

For the considered MHD problem ip is a real variable 
with asymptotic 



COs(i4' z T — (f>i), for r 



-oo, 



Df cos(K z t + 0f), for r — > +oo. 



Solving analogous 2x2 matrix problem we derive 
C^Dg = ^, C^D U = ** 



(F31) 



(F32) 



and for the phase and amplification of the signal we have 



where 



tan<pf = 

SigC u "fi Sj u Cg 

Jv (s[gS u T s^ u Sg) + (sigC u T Si u Cg) 



s ig = sin(0j - 4> g ), s iu = sin((^i - U ), 
s g = shi(0g), s u = sin(> u ), 
c,, = cosOA e u = cos(4> u ). 



(F33) 

(F34) 
(F35) 



(F36) 
(F37) 
(F38) 



- g - cos(^ g ), c 
For large enough wave-vectors we have the asymptotic 

cj>g(K 2 z » 1) = 0, MK 2 Z » 1) = -\. (F39) 



5. Random phase approximation statistical 
problem 



Finally simple angle averaging 



* — A/"(0i) = 2 - s 2 „ = 2 - V (F40) 
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FIG. 9: Numerical solution of Schrodinger equation for am- 
plitude of Alfven waves. Property function is approximation 
by two sinusoid connected with cusp corresponding of 5-like 
potential in relevant quantum mechanical problem. Initial 
conditions are: K z = 0.05 = —1000) = cos(_R' z r — tpo) ; 
d T i/>(r = —1000) = —K z svsx{K z t — tpo), and <po = § • 
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FIG. 10: Numerical solution of effective Schrodinger equation 
illustrated amplification of amplitude for small values of wave- 
vector. Initial conditions are: K z — 0.05 tp(T = —1000) = 
cos(K z t — tpo) ; d r y>(r = —1000) = — K z sva{K z T — 950), and 
tpo = TV. 



gives 



A{Kl) = --l = --l. 



(F41) 



In such a way we analyzed the relation between quan- 
tum mechanical treatment and MHD one for the effective 
Schrodinger equation. 



a. Long-wavelength approximation 

For small wave-vectors K\ <C 1 we have 5- 
potential approximation in the effective Schrodinger 



FIG. 11: Phase analysis of wave equation:the phases of odd 
(f> u and even </> g function of Schrodinger equation. Phase dif- 
ference Q gu = (j> g — ip u parameterized reflection coefficient of 
quantum mechanical problem and amplification coefficient of 
MHD problem. 




FIG. 12: Continuous curve above is phase of odd wave func- 
tion </> u , dots curve under is phase of even wave function <f> s 
shifted by f . 



equation Eq. (|F71 
1 



(1 + T 2)S 



XS(t), A 



(1 



The substitution of A — * 2A gives exacts result at long- 
wavelength limit. The wave function is continuous at 
t = i.e. 



^(-o) = v(+o), 



(F43) 



but the first derivative have a jump which can be calcu- 
lated integrating Eq. (|F7[) in a small vicinity of r = 

(F44) 



d T ?/>(+0) - d T ^(-0) = \ip(0), 



we use obvious alleviation of the notations. For the wave 
function 

ip(r < 0) = cos(K x T—<f>i), ip(r > 0) = Df cos(K z T+(f>f) 

(F45) 
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FIG. 13: Amplification of Alfven waves keystone of the theory 
of heating of accretion disks; amplification in Bell as function 
of logarithm of the dimensionless wave- vector; 1 Bell means 10 
times energy amplification. The dotted line is the 5-like long 
wavelength approximation which is adequate when the am- 
plification is significant A S> 1. This exact solution at K y = 
is a test for Monte Carlo calculations in the general case. 



one can easily solve the equations Eq. (|F43|) and 
Eq. (|F44j) which gives 



The averaging of the amplification coefficient with re- 
spect of the initial phase gives 



A = 



D 



<A<k 



A 2 



1 



2KI %K 2 Z 8 \k z V A 



ttA 



= I + — * — = - ( — — I » 1. 

(F48) 

In other words the Alfven waves amplification can be 
significant for the long waves. For the eigenfunctions in 
this approximation we have 



l 

7g 
2K Z 
A 



sin(K z |r|- 7g ), K z >0, (F49) 



4X 2 



< 1, 



7T AK 7 



2 7s ~ 2 



■ sxd.(K z t), 



7T 

2 



7g ~ 7T, 



4K« 



(F50) 
(F51) 
(F52) 
(F53) 
(F54) 



D { = J 1 + cos(20) + cos 2 (0) , (F46) 



, cos qj\ . ,„ 
6 f = arccos( — — ). (F47) 
Df 



The substitution of this approximate formula for s gu in 
the general formula for the amplification Eq. (|F41|) re- 
produces the long wavelength result Eq. (|F48[) . 



